Numerical solution of isospectral flows

In this paper we are concerned with the problem of solving numerically isospectral flows. These flows are characterized by the differential equation L′ = [B(L), L], L(0) = L0, where L0 is a d × d symmetric matrix, B(L) is a skew-symmetric matrix function of L and [B,L] is the Lie bracket operator. We show that standard Runge–Kutta schemes fail in recovering the main qualitative feature of these flows, that is isospectrality, since they cannot recover arbitrary cubic conservation laws. This failure motivates us to introduce an alternative approach and establish a framework for generation of isospectral methods of arbitrarily high order. 1. Background and notation 1.